Anomalous Hall effect in 2D Dirac band: link between Kubo-Streda formula and 

semiclassical Boltzmann equation approach 
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The anomalous Hall effect (AHE) is a consequence of spin-orbit coupling in a ferromagnetic 
metal and related primarily to density-matrix response to an electric field that is off-diagonal in 
band index. For this reason disorder contributions to the AHE are difficult to treat systematically 
using a semi-classical Boltzmann equation approach, even when weak localization corrections are 
disregarded. In this article we explicitly demonstrate the equivalence of an appropriately modified 
semiclassical transport theory which includes anomalous velocity and side jump contributions and 
microscopic Kubo-Streda perturbation theory, with particular unconventional contributions in the 
semiclassical theory identified with particular Feynman diagrams when calculations are carried out 
in a band-eigenstate representation. The equivalence we establish is verified by explicit calculations 
for the case of the two-dimensional (2D) Dirac model Hamiltonian relevant to graphene. 
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I. INTRODUCTION 



In spite of its long history as a basic ferromagnetic 
metal characterization tool, the theory of the anomalous 
Hall effect (AHE) continues to be a subject of confusion 
and debated Theoretical descriptions of the dc AHE in- 
variably involve long and complex calculations and usu- 
ally do not lead to simple general results with transparent 
interpretations^^ Theories of the AHE normally fo- 
cus on particular simple model Hamiltonians, and ignore 
interactions apart from mean-field exchange potentials 
that encode magnetic order. Even with these simplifica- 
tions, the AHE problem tends to be difficult because the 
effect usually follows mainly from density-matrix linear 
response that is off-diagonal in Bloch state band indices, 
i.e. from induced interband coherence rather than sim- 
ply from changes in Bloch state occupations. Farraginous 
theoretical results have followed from the application of 
different methods to the same models£&£& 9 Ri gorous 
approaches based on Green's function techniques, like the 
Keldysh or Kubo-Streda formalisms, have the advantage 
of being systematic but can be technically more difficult 
to implement and are often not physically transparent. 

Renewed interest in reaching consensus on a general 
theory of the AHE has been fueled in part by the real- 
ization of diluted magnetic semiconductor (DMS) ferro- 
magnets which have strong spin-orbit interactions, large 
anomalous Hall effects, and variable carrier densities 
and magnetic properties. DMS ferromagnets can serve 
as a testing ground for the different theoretical pre- 
dictions concerning the anomalous Hall effect £L±L12iM 
Interest has also been increased by the (phenomeno- 
logical) demonstration that a relatively simple intrin- 
sic contribution that is a momentum-space Berry phase 
band-structure property dominates the AHE in many 
ferromagnets^^^^^ 2 ^ 2 ^. 



From one point of view, the main reason for the phys- 
ical opaqueness of formal microscopic approaches to the 
AHE is that the Hall current is much weaker than the 
current parallel to the electric field. No Hall contribution 
appears in usual theories of the dc conductivity which 
use Gaussian disorder potential distribution models and 
evaluate conductivity to the leading order of the small pa- 
rameter £ = l/(kpl sc ), where l sc is the disorder scattering 
length. Higher order terms must therefore be considered 
- and these are unfortunately plentiful. Bookkeeping be- 
comes a challenge. Additionally typical metallic or semi- 
conducting ferromagnets have many occupied bands and 
anisotropic band structures, making it difficult in general 
to obtain analytical solutions to the relevant transport 
equations. 

The semiclassical description of transport coefficients 
through a Boltzmann equation does lead to a Hall con- 
tribution if skew scatterin g 23 ' 24 ^ 5 is accounted for in 
the collision term. There are however, other contribu- 
tions to the Hall effect that are ignored in such a the- 
ory including anomalous velocit y 12 ' 17 ! 18 ' 26 ! 27 and the side 
jumr^S&^^i^Li 3 ^ 3 . effects. The presence of these addi- 
tional contributions is directly linked to induced nonzero 
interband elements of the density matrix, in other words, 
the interband coherence, as recognized since the work of 
Kohn and Luttmger^i 3 - 4 - 

In a previous paper several of us demonstrated^ 3 , that 
by utilizing the Fermi golden rule and a newly derived 
gauge invariant expression for the coordinate shift due 
to a scattering event one can construct a semiclassical, 
physically transparent description of all contributions to 
the anomalous Hall effect and obtain the same final re- 
sults as derived in a more complicated way using the mi- 
croscopic perturbative approach of Luttingcr^ The ma- 
jor drawback of this semiclassical approach is the fact 
that the identification of contributions to the AHE is not 
systematic; it is not obvious for example that other con- 
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tributions of the same size are not being neglected and 
therefore, the assertion that we do include all main con- 
tributions must be verified by comparison to other ap- 
proaches since the semiclassical approach lacks the rigor 
of the microscopic approaches. The main goal of the 
present paper is to provide a more detailed explanation 
of the semiclassical approach and to link it with a more 
systematic and formally rigorous approach based on the 
Kubo-Streda formula. The various different contribu- 
tions that have been identified in the semiclassical ap- 
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proach, can, for the particular case of two-dimensional 
Dirac bands, all be directly linked to particular Feynman- 
diagram subsums in the Kubo-Streda formalism. The 
correspondence is summarized in Fig. [TJ in which the 
diagrammatic expansion of the Kubo-Streda formalism 
in the chiral (band-eigenstate) basis is summarized. In 
Fig. Q] Feynman diagrams are grouped in the three ma- 
jor contributions that arise naturally in the semiclassical 
approach: intrinsic, side-jump, and skew-scattering. 



intrinsic contribution 




side jump 



+ 



-T+ ^ ^^^ T" 

skew scattering 



FIG. 1: Graphical representation of the AHE conductivity within the chiral (band eigenstate) representation. The two bands of 
the two-dimensional Dirac model are labeled " ±" . The subsets of diagrams that correspond to specific terms in the semiclassical 
Boltzmann formalism are indicated. A detailed discussion of each set of diagrams is given in Sec. IIVI 
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We organize the rest of the article as follows. In SccHTl 
we introduce the model Hamiltonian for which we demon- 
strate the equivalence of the two approaches. In Sec. IIIII 
we provide a formal description of both techniques and 
discuss the approximations involved in their application. 
In Sec. IIVI we discuss different contributions identified in 
the Boltzmann approach: intrinsic, side-jump, and skew 
scattering contributions, and show their equivalence to 
partial sums of diagrams in the Kubo-Streda formalism. 
In Sec. |V]we discuss some diagrams that appear in the 
Kubo expansion beyond our approximations. Finally we 
summarize our results in Sec. IVII 



II. MODEL HAMILTONIAN 

Although most of the discussion in the following sec- 
tions is general and applicable to electronic structure in 
any number of dimensions, we choose as a concrete exam- 
ple the massive 2D Dirac model Hamiltonian with ran- 
domly placed weak <5-function-like spin-independent im- 
purities. We have chosen this model due to its simplic- 
ity which permits the derivation of compact analytical 



expressions that facilitate comparison between the semi- 
classical Boltzmann approach and the microscopic Kubo 
approach. The impurity free 2D Dirac Hamiltonian is 



H = v(k x a x + kyOy) + Aer 2 



(1) 



where a x and a y are Pauli matrices. (In the case of 
graphene, this degree-of-freedom represents the two sites 
per cell on a honeycomb lattice.) Here and throughout 
the text we take K = 1. This Hamiltonian breaks time- 
reversal symmetry and therefore has a non-zero Hall con- 
ductivity. Kubo formula calculations of the dc-limit of 
the Hall conductivity for this Hamiltonian have been al- 
ready performed 35 in the self-consistent non-crossing ap- 
proximation with the following result which can be ap- 
plied to the charge and spin Hall conductivities in metal- 
lic graphene: 
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3(vk F y i 
"(4A 2 + (Jfc F ) 2 ) 2 J 



4(vk F ) 2 
4A 2 + (i>fc F ) a 

(vk F ) 4 A 



27miV 4 (4A 2 4>fc F ) 2 ) 2 ' ( 2 ) 

where kp is the Fermi momentum and the parameters Vg 
and V\ characterizing the disorder distribution as defined 
below. 
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The non-zero mass A opens a gap in the spectrum and 
splits the Dirac band into sub-bands above and below the 
gap with dispersions 



4 =±v/A 2 + (wfc) 2 



(3) 



where k — |k| and the labels ± distinguish bands with 
positive and negative energies. In what follows, we will 
assume that the Fermi energy is positive and the Fermi 
level lies above the gap, as illustrated in Fig. O 




FIG. 2: Energy dispersion for the massive 2D Dirac Hamilto- 
nian. Filled eigenstates are represented by the colored region. 
The Fermi level is assumed to be above the gap. 

The chiral basis states that diagonalize the Hamilto- 
nian in Eq. (QJ are ip£ = (l/v / V F )e lk r |u^ : ), where 



+ _ , cos(0/2) 
1 k; ~ \ Bin(0/2)e* )' 1 k ' 



sin(6>/2) 
-cos(0/2)e^ J ' 



(4) 

cos(0) = A/y/(vk) 2 + A 2 , and tan(</>) = k y /k x . In the 
eigenstate representation (Eq.(|4|)), the disorder and ve- 
locity operators have the following matrix form 



sin(0) cos(0) 
i sin(</>) — cos(#) cos(</>) 



sin(0) sin(</>) 
-zcos((/>) — cos(#) sin(</>) 



-i sin(0) — cos(#) cos((f>) 
— sin(0) cos(4>) 

(5) 

i cos(0) — cos(9) sin(0) 



sm(6) sin(0) 



and 



(6) 



v k , k = v£ k 



(u+\u+) (u+\u k ) . _ x 

K>\<) K>K) >- Vk ' kX 



l-sin 2 (|)(l-e^-^) sin(6»)i^4 
sin(fl) 1 e 



l-cos 2 (f)(l-e^-^) , 
(7) 



We consider the model of randomly located ^-function 
scatterers: V(r) = Y^i ViS(r — Rj) with i?j random and 
strength distributions satisfying {Vi)dis = 0, (Vf 
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y 2 ^ and (V^% s = F] 3 ^ 0. This model is dif- 
ferent from the standard white noise disorder model 
in which only the second order cumulant is nonzero; 
(l^k'kl )dia = n i^o 2 where n% is the impurity concentra- 
tion and other correlators are either zero or related to this 
correlator by Wick's theorem. The deviation from white 
noise in our disorder model is quantified by V\ ^ 0, and 
is necessary to capture part of the skew scattering contri- 
bution to the anomalous Hall effect. Note however that 
within this model, as we will show in later sections, there 
are other contributions to the skew scattering that arise 
from the Gaussian disorder distribution which turn out 
to be independent of disorder strength, e.g. see Sec. lIV Cl 
In what follows we use +/— for the upper/lower chi- 
ral band indices. For example, from Eq. ([5]) = 

|2\ _ 
dis — 



oK4>-4>') 



'k'k 



uL-isin(0)-cos(0)cos(<£)] or from Eq. © (\V^ +l2 
n^ 2 |l-sin 2 (f)(l 
state basis, we will refer to the velocity matrix elements 



Working in the eigen- 
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and v i as the diagonal velocities and to v , and 

x/y ° x Iv 



v xty as t ne on? diagonal or interband velocities. 



/y 



where orbital disorder matrix element. 



III. THE BOLTZMANN EQUATION AND THE 
KUBO-STREDA FORMULA 

A. The semiclassical Boltzmann equation. 

There is a substantial literature ^ 12 ' 25 ' 33 on the appli- 
cation of Boltzmann equation concepts to AHE theory. 
However, stress was usually placed only on one of the 
many possible mechanisms for a Hall current. In this sec- 
tion we briefly review this approach, incorporating some 
new insights from recent worki^ 3 - 

The semi-classical Boltzmann equation (SBE) de- 
scribes the evolution of the electron distribution function 
as though electrons were classical particles labeled by a 
band index and living in the crystal's momentum space. 
From a quantum mechanical point of view, it is clear that 
this approach is not universally applicable because it vi- 
olates the Heisenberg uncertainty principle by having a 
momentum-distribution function at each point in space. 
A rigorous treatment of the electronic state evaluation in 
,the presence of disorder should generally consider the en- 
tire density matrix, including off-diagonal elements. For- 
mal justifications of the Boltzmann description are usu- 
ally made in terms of the properties of wave packets that 
have well defined average momentum, center of mass co- 
ordinate, spin, intrinsic angular momentum, etc£^ The 
distribution functions then acquire the meaning of the 
probability distribution of wave packets that behave in 
many respects like classical particles. 37 

The semi-classical distribution function evolves both 
due to hydrodynamic particle fluxes and due to collisions 
with impurities. Quantum mechanically, the elementary 
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scattering process is described by the scattering matrix 
from one state to another in the presence of a perturb- 
ing potential. In the semi-classical approach, scattering 
can be accompanied by changes of coordinate. As we 
discuss below, these have to be taken into account when 
constructing a kinetic equation for the semi-classical dis- 
tribution function that captures the AHE. 



Since cojy is symmetric, the leading contribution to 

appears at order V 3 . The Ujf, contribution to the scat- 
tering rate low contains symmetric and antisymmetric 
parts. The symmetric part is not essential since it only 
renormalizes the second order result of Eq. (fT4")l . The 
sum in parentheses in Eq. (| 1 5[) can be rewritten in the 
form 



1. The golden rule. 

The golden rule connects the classical and quantum 
descriptions of a scattering events It shows how the 
scattering rate loih between states with different quan- 
tum numbers I and I' is related to the so called T-matrix 
elements: 

uvi =2Tr\T l , l \ 2 5(e l , - e,). (8) 

The scattering T-matrix is defined as 

T n = {l'\VHi) (9) 

where V is the impurity potential operator and is 
the eigenstate of the full Hamiltonian H = Ho + V that 
satisfies the Lippman-Schwinger 



IVO = 10 + 



v 



ei - Hq + ir\ 



(10) 



For weak disorder one can approximate the scattering 
state |V>z) by a truncated series in powers of Vw = {l\V\V) 
and take I — (k, fi) as the combined (momentum, band) 
index of the eigenstate \l) = |k, /z) of the unperturbed 
Hamiltonian H : 



T 

i" 



-\n + 



(ii) 



For example the T-matrix up to the second order in V is 



■jf ei - ei" + it] 



(12) 



Substituting Eq. (|12p into Eq. (jHJ) one can find the scat- 
tering rate up to the 3rd order in the disorder strength 



where 



(2) , (3) , 



4?) ^27T(\V w \ 2 ) dls S(e l -e l ,) 



(13) 



(14) 



(3) o I {Vll'Vl'l"Vl"l)dis . \ x , \ ne\ 

lo,,, — Z7r | > ; h c.c. \ 6{ei — til). (15) 
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The skew scattering contribution to the Hall effect fol- 
lows from the asymmetric part of the scattering rate^ 



(a) — urn 

Will = « • 

u 2 



(16) 



2R c (V ll ,V l , lll V ll , l ) dl 
ej-ej// 



(17) 



27r£,„<5(e i - e V i)lm{VwVii V iVii 



l J dis 



The first term is symmetric under the exchange of indexes 
I «-> V (note that this needs the fact that ei = e;/ due 
to the delta- function in Eq. (fl~5|) ). hence this sum does 
not contribute to the asymmetric part of the scattering 
amplitude and one can concentrate on the second one: 

Lo ( * a) = ~(2TifY f 8{e l - e,n)hn(V w V vvl V m ) dia S(e l - e v ), 

with the superscript 3a meaning that this is the anti- 
symmetric part of the scattering rate calculated at order 
V 3 . 



2. The coordinate shift 

In a semi-classical description of wavepacket motion in 
a crystal, scattering produces both a change in the direc- 
tion of motion and a separate coordinate shift £ ' 28 i 29 ' 30 i 33 
In the lowest order Born approximation for a scalar dis- 
order potential one can derive an expression for the co- 
ordinate shift which accompanies scattering between two 
band states. The shift does not depend explicitly on the 
type of impurity and can be expressed in terms of initial 
and final states onlyi^ 3 - 3 - 

d d 
5v V i = \uii\i—up) - (ui\ij^ui) - D k / )k arg[(u//|ttj)]. 



dk 



(19) 

where arg[a] is the phase of the complex number a and 



D 



k '. k 



d d_ 
dk> + dk 



The topological interpretation of Eq. (jT9|) has been ex- 
plained in Ref . [H. The first two terms in Eq. (fT9"|) have 
been known for a long time. The last term was first 
derived only recently in Ref. [H but is an essential con- 
tribution which makes the expression for the coordinate 
shift gauge invariant. 



Kinetic equation of the semi- classical Boltzmann 
distribution 



Eqs. |(8J) and (|19|) contain the quantum mechanical in- 
formation necessary to write down a semi-classical Boltz- 
mann equation that takes into account both the change 
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of momentum and the coordinate shift during scatter- 
ing in a homogeneous crystal in the presence of a driving 
electric field E. Keeping only terms up to the linear order 
in the electric field the Boltzmann equation reads^ 3 . 



de. 



= - Ez< [fl 



-fv 



where Voz is the usual velocity 

vqi = dei/dk. 



(20) 



(21) 



Note that if only elastic scatterings with static impuri- 
ties are responsible for the collision term, the rhs. of 
Eq. (|20|) is linear in /; and not in fi(l — /;/). This 
property follows from the fact that Heizenberg time evo- 
lution equations for creation or annihilation operators in 
a noninteracting electron system are linear and can be 
mapped to the Schrodinger equation for amplitudes in 
a single particle system. (For further discussion of this 
point see Appendix B in Ref. [3^ .1 

This Boltzmann equation has the standard form except 
for the coordinate shift effect which is taken into account 
in the last term in the collision integral on the rhs of Eq. 
(|20| . This term appears because the kinetic energy of 
a particle in the state /' before scattering into the state 
I is smaller than e; by the amount eE • Srw = — eE • 
Sri' i. The collision term does not vanish in the presence 
of an electric field when the occupation probabilities fi 
are replaced by their thermal equilibrium values because 



fo(ei) - fo(ei - eE • 6r w ) ~ - 



eE • St vi ^ 0. The 



last term in the collision integral in Eq. ([20| 6 accounts 
for this interplay between coordinate shifts and spatial 
variation of local chemical potential in the presence of an 
electric field. 

The total distribution function /; in the steady state 
(dfi/dt — 0) can be written as the sum of the equilibrium 
distribution /o(e;) and non-equilibrium contributions, 



fl = h{ei)+9i+gT 



(22) 



where we have split the non-equilibrium contribution into 
two terms gi and g adzst which solve independent self- 
consistent time-independent equations) 3 ^ 



eE ■ v f 



dfojei) 
de. 



- {gi -gv) 
i 



(23) 



and 



V 



vu> gi 



^adist 



9i 



adist 



-a/ofo) 

de. 



eE • 8t v 



(24) 

(The label adist stands for anomalous distribution.) A 
standard approach to solving these equations in 2D is to 
look for the solution of Eq. (|23|) in the form (see e.g. Ref. 
22 



9l = (" 



dfojei) 
dei 



eE • (A^vqi + B^v i x z) 



(25) 



where z is the unit vector in the out-of-plane direction. 
We will assume that the transverse conductivity is much 
smaller than the longitudinal one so that >> B^. 
One then finds by the direct substitution of Eq. 
into Eq. ([20]) that 



-4, 



where 



(26) 



(27) 



and <{> and (/>' are the angles between vo; or Voz' and the 
x-axes. This completes the solution of Eq. (|23[) for gi. 
We will show in Sec.III.B.l how a similar ansatz can also 
solve Eq. (|24|) for the anomalous distribution g® dlst . 



4- Modified velocities and the anomalous Hall effect. 

To find the current induced by an electric field and 
hence the conductivity we need to derive an appropri- 
ate expression for the velocity of semiclassical particles, 
in addition to solving the SBE for the state occupation 
probabilities. In considering the AHE, in addition to the 
band state group velocity v / = dei/dk one should also 
take into account velocity renormalizations due to the 
accumulations of coordinate shifts after many scattering 
events (the side-jump effect) and due to band mixing by 
the electric field (the anomalous velocity effect) ; 30 i 33 



dk 



(28) 



The second term in Eq. (1281) captures changes in the 
speed at which a wave packet moves between scattering 
events under the influence of the external electric field 39 
and Fi is the Berry curvature of the band2& 



(Fi) k = eijklm 



, dui dui 



(29) 



The last term in Eq. (|28[) is due to the accumulation of 
coordinate shifts after many scatterings. Combining Eqs. 
(l22l) and (l28l) we obtain the total current 



/V/ 



(30) 



5. Mechanisms of the AHE. 

From (|3"0"]) we can extract the off-diagonal Hall con- 
ductivity which is naturally written as the sum of four 
contributions: 



total rr^ n ^ 

xy xy 



adist 
xy 



xy 



sk 
xy 



(31) 
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The first term is the so called intrinsic contribution 



■ XI) 



(32) 



which is due to the anomalous velocity of free electrons 
under the action of the electric field 



Ft x eE. 



(33) 



This term is called the intrinsic contribution because it 
is not related to impurity scattering, i.e. its evaluation 
does not require knowledge of the disorder present in the 
system. The intrinsic AHE is completely determined by 
the topology of the Bloch band. 

The next two contributions follow from coordinate 
shifts during scattering events: 



adist 
7 xy 



(34) 



is the conductivity due to the anomalous correction to 
the distribution function while the direct side-jump con- 
tribution 



(35) 



is due to the side-jump velocity i.e. due to the accumu- 
lation of coordinate shifts after many scattering events. 
Since coordinate shifts are responsible both for <J X y lst and 
for a^y there is, unsurprisingly, an intimate relationship 
between those two contributions. We will demonstrate 
their equality by evaluating a X y lst below. Because of this 
relationship, in most of the literature, cr%y lst is usually 
considered to be part of the side-jump contribution, i.e. 
(Txy lst + °xy ~ * a xy Here we will distinguish between the 
two when connecting the SBE approach to the Kubo- 
Streda formalism. 

Finally, a s x k y is the skew scattering contribution 



e22{gi/Ey){v a i)x- 
i 



(36) 



There is potentially one more contribution which has 
not been included in Eq. (|3ip and which involves the 
product of g adlst and the side jump velocity w/'/<5r//;. 
Such a term does not contribute to the Hall conductiv- 
ity in the present model calculation of isotropic scatter- 
ing and band structure. It may however, in principle, 
give rise to a non-zero contribution given an appropriate 
anisotropic scattering and anisotropic Fermi distribution. 
We do not consider this term further here. 

In Sec. IV we evaluate each of these four contributions 
for the particular example of the 2D Dirac model. 



B. The Kubo-Streda formula. 

A fully quantum mechanical formally exact expression 
for the conductivity is provided by the linear response 
Kubo theory. There are various equivalent formulations 
of this approach. In one version the conductivity at zero 
temperature T is expressed in terms of causal Green's 
functions. This formalism has previously been applied to 
the AHE in an electron system with Rashba spin-orbit 
coupling. 5 In the present paper we work mainly with an- 
other version of the Kubo formula, derived specially for 
the the dc conductivity. Its advantage is that one should 
not perform the to — > limit at the end of calculations. 

The starting point in this approach is the Bastin 
formula^ 



e 2 f°° de 
-»V / —.f{e)TT[-vJ{e-H)v 
V 2m 



dG R 
J ~dT 



(39) 



which can be manipulated to give the Kubo-Streda for- 
mula for the dc T = Hall conductivity a xy = a x 



1(a) 



1(b) 
>xy 



all where: 



li a) = ^TV^G^KG^)^ (40) 



The skew scattering contribution is independent of the 
coordinate shift and of the anomalous velocity. It is non- 
zero for an isotropic energy dispersion e/ only when the 
scattering rate bJw is asymmetric. Writing the trans- 
verse velocity in the form (vqi) x = \voi\ cos(0), then using 
Eq. |J57J| and substituting Eq. (QSJ) for g t into Eq. 
one finds 



< = eEJ(f^(^/^)l«oi|cos(^) 



At zero temperature this expression simplifies to 



■xy 



e 2 {rl r ) 2 VF^k Fl , 



(37) 



(38) 



yl(b) 



II 



AttV 



Tr(v x G R (e F )v y G R (e F )+v t G A (e F )v 3 G A (e F )) c 



(41) 



4irV 



def(e)Ti[v x G R (e)v. 



G R (e) 



GR{e) -v y G R (e)-v x G A (e)vy 



" de 
G A {e) 



de 
G A (e) 
de 



dc 



v y G A (e)} = ec 



dN(e 



dB 



(42) 



e=e F ,S=0 



here the subscript c indicates a disorder configuration 
average. The last contribution, a xy , was first derived by 
Streda4^ In these equations B is an external magnetic 
field that couples only to the orbital degrees of freedom, 



7 



and N(E) is the integrated density of states of the dis- 
ordered system and G R ^ A (e F ) — (e F — H ±i5)~ 1 are the 
retarded and advanced Green's functions evaluated at the 
Fermi energy. Every term in a^. y depends on products of 
retarded Green's functions only or products of advanced 
Green's functions only. Since every such term has poles 
on the same side of the imaginary plane, disorder correc- 
tions to this contribution are small because they contain 
at least one (small) factor niV 2 which is not compensated 
by energy and momentum integrations. Hence only the 
disorder free part of a^ y ^ s important in the weak dis- 
order limit, i.e. this contribution is zeroth order in the 
parameter £ = l/k F l sc . It can also be shown by a similar 
argument that in general cr^, is of order £ and can be 
neglected in the weak scattering limit. 4 - Thus, important 
disorder effects are contained only in cr£° . In what follows 
we neglect a™ and retain only the clean limit ajy con- 
tribution. The following paragraphs concern only a 1 ^. 

The effect of disorder on the disorder-configuration av- 
eraged Green's function is conveniently captured by the 
use of the T-matrix, defined by the integral equation 
T = W + WG T where W = V S{r - n). From 
this one obtains 



result for the 2D Dirac model Hamiltonian is 

„ (1 + cos 2 6) ( sin 2 6 \ 

^ = 8VT C ° S 9 l(lTwf ffl+ V + > + 3cos 2 0) ) a 

(47) 

where T = 7m 2 V^ 2 /4v 2 . Incorporating this result in cr^ 
we obtain for e F > A 

a xy = Ek Tr [voxGvyG] 

(48) 

_ e 2 cos 8 ( i , 3 sin 2 , 4 sin 2 8 \ 
in y 1 (1+3 cos^ 8) "r (l+3cos a 0) 2 J 

This is equivalent to the part of Eq. ([2]) that is inde- 
pendent of V\. In order to understand more fully the 
relationship between this result and the more physically 
transparent SBE approach, we need to analyze the Kubo 
approach in the chiral band eigenstate basis. As de- 
scribed in the remaining part of this section, we are able 
to make the connection very clear. 

We have also obtained the same result using the time 
ordered Green's functions formalism. The equivalence of 
the two results is shown in detail in Appendix [B] 



G = Go + G TGo = Go + G ZG. (43) 
Upon disorder averaging we obtain 

£ = (W) c + (WG W) C + (WG a WG Q W) c + ... (44) 
To linear order in m this translates to 

E(z,k) = n i y k , k + ^2^ lk ^o(k' > «)V k / jkj (45) 

k 

with Hc.k' = V(k — k') being the Fourier transform of 
the single impurity potential, which in the case of delta 
scatterers is simply Vq. 

Within the T-matrix formalism in the normal spin- 
basis we will take V\ = for simplicity since, as it will be 
shown when considering the Kubo formalism within the 
chiral basis, only a couple of diagrams contribute to part 
of the skew scattering for V% ^ and can therefore its 
contribution can be computed separately (see Sec. lIVC]) . 

Replacing the bare Green's function by the self con- 
sistent Born approximation Green's function, Eq. (|43[) . 
in the expression for a£y does not capture all the con- 
tributions to lowest order in £. One must also compute 
the ladder diagram correction to the bare velocity vertex 
defined by V a {ep) = v a +5v a (eF), where 5v a (eF) satisfies 

Sv a (e F ) ^^G fl ( ff KG A M 

k 

+ ^_Y^G R {e F )8v a {e F )G A (e F ). (46) 

k 

The explicit calculation of this vertex correction in the 
<j z eigenstate basis is described in Appendix [X] The final 



1. Kubo-Streda approach in the band eigenstate basis. 

In the chiral Bloch eigenstate basis, specified in Eq. 
0, the disorder free Green's functions in momentum 
space can be written in the form 

G n/A {e) = K){uj\ + 

e — eZ i ifj e — e k ± irj 

From a calculational point of view the bare Green's func- 
tion is simpler in this basis, but we must pay the price in 
the disorder matrix of the Hamiltonian which will con- 
tain, as seen in Sec. II, intraband (V^T, Vj^r) and in- 
terband (V^J , V^J~) matrix elements. 

In the weak disorder limit, the energy width of Bloch 
state spectral peaks is smaller than the gap, allowing 
us to ignore direct interband scattering. This does not 
mean, however, that interband matrix elements of the 
disorder potential are not important and can be disre- 
garded. In the rest of the text we will show that interband 
disorder matrix elements, and the interband coherence 
they imply, are in fact crucial to fully explain disorder 
contributions to the AHE. Although interband scatter- 
ing probability remains zero, interband disorder matrix 
elements produce virtual transitions that mix states in 
the two bands in a way which is ultimately relevant. 

When evaluating diagrams in the chiral basis, disorder 
lines associated with interband and intraband scattering 
processes can be distinguished. Their separate contribu- 
tions will be easily linked to the semiclassical interpreta- 
tion of the scattering process. As we will see below, the 
relevant vertex correction infinite subsum involves only 
intraband elements V^J' and does not produce unusual 
effects in the sense that it only renormalizes the quantum 
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life time and the diagonal velocity vertex, v ++ — > T ++ . 
In contrast, vertex correction diagrams which contain 
V^7 , V^t and therefore connect Green's functions with 
different band indices, do not have to be summed to infi- 
nite order and do capture the unusual transport physics 
which leads to AHE contributions. Only a few low order 
(in terms of £) diagrams, which depend on or V^f , 
have to be included in the calculations of a xy . Note that 
these terms still contain the renormalized T ++ velocity 
vertex but the interband component does not contribute 
to this renormalization as we will see below. 




FIG. 4: Graphical representation of the vertex equation. 



3. The vertex correction. 



2. The self-energy in the non-crossing approximation. 

We consider first all possible diagrams that do not con- 
tain an interband part of the disorder potential. This 
part of the calculation is essentially the familiar standard 
calculation because the Green's functions remain diago- 
nal in band indexes and the calculations are reduced to 
single band ones. In this case the non-crossing approx- 
imation leads to a finite life time and the diagonal ve- 
locity vertex is renormalized after the summation of all 
non-crossing disorder paring lines. 

Only the imaginary part of the self-energy is impor- 
tant in the ( <C 1 limit. The self-energy diagram that 
contributes to its diagonal element for the Gaussian dis- 
order model is shown in Fig. [3l Because G R ~ (e) has a 
negligible imaginary part one can keep only the G R+ (e) 
part of the Green's function in calculations of T, R , i.e. 
Y, R = J2k> ^kk^^ fl+ ^k^k ■ Hence, the renormalized re- 
tarded Green's function is given by a standard result 

G*+ (e) - G«+ (e) = (50) 

where 

l/r+ = -2Im(S«)=2 7 r J 0-\vR\*6{e- e+ ) (51) 

Note that we have so far ignored the off-diagonal part of 
the self-energy which will contribute later when consid- 
ering the diagrammatic expansion of the conductivity in 
powers of £. It is convenient to organize the calculation 
in this way because the off-diagonal parts do not need 
to be partially summed to infinite order and only a few 
terms contribute to the desired order. 



/ \ 
/ \ 

£ R = * » k 

FIG. 3: Self-Energy diagram 



Every disorder line that connects retarded and ad- 
vanced Green's functions adds a factor njVg 2 ~ £ mul- 
tiplied by a product G R G A . While the former factor 
gives a power of £, the product of two Green's functions 
gives a £ contribution if the pair of Green's functions 
over which the momentum sum is performed are in the 
same band. Thus one can approximate^ 2 - 

G R+ {e)G A+ (e) w 27rr+^(e - e k ). (52) 

Since r + is proportional to l/£ this product of 
Green's functions compensates the small factor due 
to disorder vertexes. Therefore terms Vy + and 

Ek' Ku GR+v i +GA+v v£ are of thc samc order and 

one should perform the summation over all ladder dia- 
grams that include only Green's functions in the positive 
energy band. These ladder diagrams simply renormalize 
the diagonal matrix element of the velocity vertex oper- 
ator. 

The bare diagonal velocity is Vy + — v sin(0) sin (</>). 
The vertex equation is 

T++ = <+ + / |^+|^ k +| 2 G*+T++ (53) 

Eq. (|53p for the vertex correction is represented graphi- 
cally in Fig. Q] This equation has a solution of the form 
T+ + = av sin(#) sin(<^>) where a = r* r /r + , 

1/r*- = 2tt / $fc\V++\ 2 (l - cos(k,k'))<5( £F - e +) 

_ (t>fc F ) 2 +4A 2 
~~ 4T«((-ufc F ) 2 +A 2 ) ' 

(54) 

l/ r « = !^£ (55) 

Vf 

and vp is the Fermi velocity vp = (<9e + /dk)k=k F = 
v 2 kp/eF. Combining this result with Eq. (fB"2")) gives the 
useful weak-scattering identity 

< + <^+ - G«+T++G^+ » vtp^5{e e+). 

(56) 



4- Expansion in the number of virtual interband transitions. 



After summing the series of diagrams that contain only 
diagonal disorder vertices one can look at the effect of 
the inter-band disorder matrix elements V^7 and Vj^/ . 
Any disorder average line with this component still adds 
a small factor n.$ Vq , however now it is followed by at least 
one Green's function in the lower band G R ~ or G A ~ . The 
new Green's function products that accompany the inter- 
band disorder interaction lines, products like G R ~G A+ , 
etc., no longer can be considered as large. They are at 
most zeroth order in the parameter £ in the final expres- 
sion rather than being proportional to Since any 
part of a diagram containing either G R ~ or G A ~ adds a 
factor of the small parameter £, these diagrams do not 
have to be summed to infinite order. It is enough to 
consider only the finite set of them that contribute at 
order £°, the set of diagrams with only one G R ~ or G A ~ 
line. The full set of such diagrams in the non-crossing 
approximation is shown in Fig. Q] and in individual sets 
in Figs. [H El andfTTl Of course, to find the total con- 
ductivity, one should consider them together with any 
number of intra-band disorder averaging lines which are 
already included by replacing v ++ by T ++ . 



IV. IDENTIFICATION OF SBE AND KUBO 
FORMALISM AHE CONTRIBUTIONS. 





FIG. 5: Diagrams, without disorder paring lines, contributing 
to the Hall current. Plus/minus refer to the band index with 
positive/negative energy 



Hamiltonian we find 

F ± (k) = 2Im(d ky u±\d kx u±) = T 



Av 2 



2(A 2 + (vk) 2 ) 3 / 2 



(57) 

where ± in F ± (k) refers to upper/lower bands. Since the 
intrinsic contribution comes from the whole Fermi sea, 
both the totally filled negative energy and the partially 
filled positive energy bands contribute to it. Substituting 
Eq. (JSU) into Eq. J32J we find 



intr _ 2 f°° d a k p-(]A 
"xy — e Jo l27p r \ K ) 

_ e 2 r 4*F+{k) = j£± 

Jk<k F (2tt) 2 \ > 4n^/A 2 + (vk F ) 2 

which coincides with the first term in Eq. ([2]). 



(58) 



After discussing the two formalisms we are now ready 
to demonstrate their equivalence and to identify the set 
of chiral basis diagrams that correspond to the various 
physical effects that have been identified in the more 
physically transparent semiclassical treatment. It is im- 
portant to note that the same result is obtained by doing 
the calculation in the normal spin basis, as shown in Ap- 
pendix El but in that case a given Feynman diagram can 
contain contributions from more than one of the effects 
identified in the semiclassical theory. Calculations in the 
band eigenstate basis are necessary to interpret the SBE 
physics because the SBE approach starts by ignoring in- 
terband coherence and adds coherence effects by hand by 
positing side-jump and anomalous velocity effects. 



A. The intrinsic contribution 

1. The intrinsic contribution from the SBE 

According to Eq. (p?2"|) , to find the part of the Hall con- 
ductivity not related to impurity scattering one should 
first find the Berry curvature of both bands. In two di- 
mensions, only the z-component of the Berry curvature 
is meaningful. From Eqs. (H]) and (j2"9")l for the Dirac 



2. The intrinsic contribution from the Kubo formalism 

Since the intrinsic contribution to the AHE, as iden- 
tified in the semiclassical approach, does not depend on 
impurities, it is expected that this contribution must cor- 
respond to the Kubo formula contribution from diagrams 
that do not involve any disorder lines, i.e. the bubble di- 
agram conductivity with all Green's functions bare. In 
Fig. \5\ we show the diagrams from o l xy that are non-zero 
even without disorder. The unperturbed Green's func- 
tion is diagonal in the chiral basis given by Eq. (|3]); the 
band label in the diagram is "+" for G R+ + and "— " 
for G . The contribution to the Hall current due to 
these diagrams is 

I(int) _ e 2 r d 2 U + ~v~ + < ") r / +x _ 

a w - 2¥ J {5¥F 7^7^ °K e F- e k)- 

= e 2 A 

47T % /A 2 + (t)fc F ) 2 

(59) 

Surprisingly this result coincides with the semiclassical 
result obtained by evaluating Eq. even though we 

have not yet included the disorder free part of tr^, which 
is also non-zero in general and must be included as a part 
of the intrinsic conductivity. It turns out however that 
for the case of the 2D Dirac model all vanishes when 
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the Fermi level is above of the gap. To see that we first 
observe that a^ y is exactly the sum of two contributions, 
one of which is the same as the intrinsic conductivity 
given by the SBE Eq. while the other is equal to 

the intrinsic part of ovL but with the opposite sign: 



where F + is the Berry curvature of the upper band, de- 
fined in Eq. (|57j) . The average side-jump velocity can 
be obtained by multiplying the transition rate and the 
side-jump associated with a particular transition as in 
Eq. ([!>]): 



7 H 

'xy 



J 



F+(k)- 



(*)) 



I(int) 
Txy 



(60) 



Details of the calculation which lead to this result are 
presented in Appendix [C] This result is of course not 
surprising when we recognize that Eq. expresses the 
linear response of a clean system, i. e. the full a xy in the 
absence of disorder under the influence of an accelerating 
electric field; since a xy must equal a xy +a xy it follows that 
a ly — °xy — &i y ■ This observation also means that there 
is no special intrinsic contribution arising from the Fermi 
surface, where occupation numbers change because of ac- 
celeration, because the corresponding terms in a xy and 

all, cancel and the final result coincides with the semi- 
y 

classical prediction Eq. (|58|) . coming from all electrons 
below the Fermi surface. 

We have already calculated the sum of the two integrals 
in Eq. (j6"0|) and aiy^ (see Eqs. [3H1 EHJ) and showed that 
they have the same value. It follows that in the metallic 
regime, when ep > A, a xy is identically zero. In the case 
of the 2D Dirac Hamiltonian, the Fermi surface contri- 
bution axy int ^ = vanishes when the chemical potential 
is in the gap (—A < ep < A), but 



(0/K)^W = -£^O (61) 

Il(gap) 



(in units h — 1). Thus a xy ' is responsible for the 
quantum Hall effect in the Dirac band 
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B. Side-Jump related contributions. 

1. Semiclassical side-jump conductivities. 

Since direct inter-band scattering is not energetically 
allowed in the weak disorder limit of our model, we con- 
sider only the coordinate shift effect of scattering in the 
upper band: 



d 



d 



-D k ', k arg((u+,|u+)) 



<frk',k 

" Ok' " ' " CK " " " ' " 

(62) 
for the 

positive energy band leads to the following expression af- 
ter noting that the absolute magnitude of the momentum 
is conserved: 



where Dk'.k = -^r + Evaluation of Eq.( 



(frk'.k = 



F+{k) i x (k'-k) 

\K,K)V 



(63) 



(k) = J2k' '^k'k(* r k',k)a! 



2^o 2 / m?(-F + ){k' y - k y )S(e+ - e+) = F+k y /r«, 

,( 64 ) 

where we have used that in the lowest Born approxima- 
tion Wk'k = 271-11^ kl 2( ^( e k ~ e k')- The side-jump ve- 
locity does not produce any contribution to the total 
current in equilibrium, but in an external electric field 
the nonequilibrium correction to the distribution func- 
tion, Eq. (|25| . results in a non-zero net current. Disre- 
garding a small antisymmetric part the distribution func- 
tion correction is gk = (— dfo/de^) eE y (v^') y T tr , where 
) y = de^/dk y . It follows that the side-jump contri- 
bution to the SBE Hall conductivity is 



'xy 



e 2 A(vk F ) 2 



(65) 



2TTy/(vk F ) 2 +A 2 ((vk F ) 2 +4A 2 ) ' 



The second side-jump effect discussed in Section III 
follows from the change of energy of the scattered particle 
after the coordinate shift at a scattering event in the 
presence of an external electric field. It results in the 
anomalous correction to the distribution function that 
can be determined from Eq. (|24[) . In our case 



E. . / adist adist 

Wk,k' \ 3k - ffk' 

k' 



-dfo 
de+ 



0. 



(66) 



Substituting <5rk',k from Eq. f|63[) and looking for a so- 
lution of Eq. (|66[) in the form g^ lst = ^k x we find that 



z^eE y T tr /r q . The corresponding contribution 
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to the Hall conductivity is 

—adist „ P d 2 k ( ~ adist j rp V,,H — h 

xy - e J(2^F(5k l^vFx : - 

e 2 A(vk F ) 2 

_ 2TT^(vk F ) 2 +A 2 ((vk F ) 2 +4A 2 ) 



(67) 



For the 2D Dirac Hamiltonian we find that a®?, = ail,, 

r—i I — I y •Ly 1 

in an agreement with Ref. [6| and 1331 , where we showed 
that in the absence of direct interband scatterings these 
contributions must be identical. Thus the total side-jump 
related conductivity is just twice the one in Eq. (|67p : 



si i adist 
^ xy 1 xy 



e 2 A(vk F ) 2 



TTy/(vk F ) 2 + A 2 ((vk F ) 2 + 4A 2 ) 



(68) 



This result coincides with the second term in the total 
conductivity Eq. ^j. 
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FIG. 6: Diagrams contributing to the side-jump current 




FIG. 7: Diagrams contributing to the current that arises from 
the anomalous distribution. High symmetry with side-jump 
diagrams in Fig[6] explains the equality of two contributions 
to Hall current. 



2. Side jump effects in the Kubo formalism. 

If we look at the expression for the conductivity cr^° ~ 
Ti(v x G R v y G A ), we find two velocity operators. Tracing 
the history of how this expression was derived, the op- 
erator v x appeared when we were looking for its average 
(i> x ) and the operator v y appeared from the coupling of 
the charge current to the electric field. Thus the right 
velocity vertex in Feynman diagrams is associated with 
the immediate effect of the electric field, such as the ac- 
celeration of wave packets (due to the coupling of electric 
field to the diagonal velocity) or the mixing of states from 
different bands (due to the coupling to the off-diagonal 
velocity) . In turn, the left velocity vertex shows what av- 
erage velocity is finally induced due to this effect if this 
diagram is nonzero. 

This observation allows us to look for the diagrams re- 
sponsible for various contributions to the AHE. The side 



jump current is due to the accumulation of transverse 
coordinate shifts. These shifts are anomalous, i.e. they 
arise from inter-band coherence induced by impurities. 
In the operator language, they appear because the ex- 
pectation value of the off diagonal part of the velocity 
operator becomes nonzero on average. Hence the dia- 
grams responsible for the side jump effect are expected to 
have off diagonal elements v+~ or v~ + on their left side. 
In the SBE the side jump velocity becomes nonzero on 
average only after the electric field induces the nonequi- 
librium correction to the distribution function g^. That 
correction is due to the coupling of the electric field to 
the diagonal part of the velocity. Hence side jump di- 
agrams should have only diagonal velocity Vy + on their 
right hand side. There are four different diagrams, shown 
in Fig. [5] with such properties. Their calculation, using 
Eq. (|50f and ([56]) , leads to the same result as the one in 
Eq. The band labels at the right vertex in these 

diagrams are near the Fermi energy in the occupied band. 

The identification of the diagrams which contain the 
contribution obtained in the SBE due to the anomalous 
distribution that arises from the side-jump, <J X y ts , can 
be identified by similar reasoning. For these contribu- 
tions, the electric field is coupled to the coordinate shift 
(and hence to the off-diagonal velocity). The coupling- 
changed the band energy upon scattering, leading to the 
anomalous correction to the distribution function. After 
this the average of the diagonal velocity from such a dis- 
tribution correction becomes nonzero. Hence the corre- 
sponding diagrams should have the off diagonal velocity 
on the right and the diagonal velocity on the left. There 
are again four such diagrams in Fig. [7] whose evaluation 
leads to the same result as Eq. (|oT|) . Symmetry relations 
between the Feynman diagrams shown in Fig. [6] and in 
Fig.[7]explain why the two contributions are equal to each 
other. However, within the semiclassical formalism, the 
seemingly exact relation between the two contributions 
is not clear to us. 

We emphasize once again that here only four diagrams 
for each of the two contribution have been evaluated to 
obtain the result of order O(£ ) and not a partial infinite 
sum. (This counting assumes that T ++ has already been 
evaluated separately.) We also emphasize that although 
the explicit demonstration of equivalence between these 
specific Feynman diagrams and the side-jump calcula- 
tions which we have identified on physical grounds is for 
the Dirac model alone, the considerations on which the 
identification has been made is general, and we believe 
that the identification is generally valid. 



C. Skew scattering 

1. Skew scattering in the SBE. 

The skew scattering contribution originates microscop- 
ically from asymmetry in the scattering rate. This asym- 
metry can be estimated for known disorder scattering 
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(a) 




+ 




FIG. 8: Lowest order non-zero terms in the Born series re- 
sponsible for the asymmetric part of |T kk /| 2 . 




(b) 



(c) + 



(d) 



(e) 



(f) 





FIG. 9: Terms contributing to the antisymmetric part of 
l^k'kl 2 m the white noise non-crossing approximation at 4th 
order in the disorder potential. 



sources by truncating the Born series expansion of the 
T-matrix at low order. As discussed in Sec. IIII1 in the 
lowest Born approximation the scattering rate is sym- 
metric with respect to an exchange of indexes. The first 
non-zero contribution appears at the next (3rd) order. In 
Eq. (fl~5|) the terms in parentheses represent the absolute 
value of the squared T-matrix element |T k! k'| 2 to the or- 
der V 3 . In Fig.[8]we show the diagrams representing both 
such terms from Eq. (|15[) . This contribution vanishes for 
a Gaussian (white noise) disorder distribution. However, 
one should generally keep this contribution, even if V\ is 
small, because it leads to a parametrically different and 
hence experimentally distinguishable contribution to the 
AHE. 

At 4-th order in V the number of terms that contribute 
to the asymmetric scattering rate, and hence skew scat- 
tering, is non-zero even for Gaussian disorder models. In 
Fig. [5]we show schematically all the terms that contribute 
at 4-th order in the non-crossing white noise approxima- 
tion. For example, the diagram in Fig. [9]Ja) corresponds 
to the expression 



i 



T:0'kk"K:''kW" 



1 



e - e k „ - in 



The summation of terms in Figs. [8] and [9] leads to the 
following expression for the asymmetric part of the scat- 
tering rate: 



where 



(3a) 
k,k' 



,0) 



(3a) 



,(4a) 
J k,k' 



2[{vk) 2 



(69) 



(70) 



FIG. 10: Diagrams contributing to 3rd order skew scattering. 



2^2 



(4Q ) _ 37rAfc 2 (n t F 2 ) 



k,k' 



A[(vk) 2 + A 2 ] 3 / 2 



The dependence of w k 3 k ^ and 



(71) 



k,k' 



,(4a) 



on the impurity 

concentration n, provides valuable information on the 
physics behind asymmetric scattering. Although w k 3 k ) 
is proportional to rii and is therefore due to separate in- 
coherent contributions from all impurities, the 4th order 
Gaussian contribution w k 4 k / is proportional to nf and 
can be interpreted as due to interference between dif- 
ferent scattering centers. Given the scattering rate one 
can calculate the transverse scattering time and the Hall 
current using Eqs. ([27f and (j37)) . The skew scattering 
contribution to the Hall conductivity is then 



e 2 3A(t>fc F ) 4 



(vkf 



■xy 



47T^/(vk F ) 2 +A 2 [4A 2 +(vk F ) 2 ] 2 2imiVfi {AA 2 + (vk F ) 2 ) 2 

(72) 

This result coincides with the last two terms in Eq. ([2]). 
Note also that the second term in Eq. (|72| is inversely 
proportional to the impurity concentration rii and thus 
the transverse conductivity due to this skew scattering 
becomes large in the clean limit whenever V\ ^ 0. In 
contrast, the first term in Eq. (|72[) appears even for the 
Gaussian disorder correlations and it is independent of 
the impurity concentration. This is a result that has 
not been previously noticed: the skew scattering contri- 
bution, as defined through the collision term in the SBE, 
can give a result independent of the impurity concentra- 
tion and disorder strength. Of course this contribution is 
subdominant if disorder is weak. 



2. Skew scattering in the Kubo formalism. 

We have shown above that intrinsic and side-jump re- 
lated diagrams contain at least one off-diagonal velocity 
operator element, v^. Skew scattering is, from this con- 
ceptual point of view, the most conventional contribution 
to the AHE since it depends only on velocity operator 
matrix elements that are diagonal in the band indices, 
v x/yi 01 more accurately its renormalized version Tj^. 
In Figs . [TU1 and [TT1 we show respectively the leading order 
diagrams for non-Gaussian correlations and for Gaussian 
disorder in the non-crossing approximation. 



13 




FIG. 11: Leading diagrams in self-consistent non-crossing ap- 
proximation for the skew scattering contribution to the AHE 
in the semiclassical approach. 




FIG. 12: A diagram with crossing disorder lines. 

The correspondence between calculations of skew scat- 
tering contributions based on the Kubo formalism and 
the SBE becomes obvious when we compare Figs. Upland 
[11] with terms contributing to the asymmetric scattering 
rate in Figs. [8] and [9] Comparing both sets of figures 
shows that the central part in Figs. ITOl and ITT1 coincides 
with the formal perturbation expansion of the asymmet- 
ric part of |Tk.k'| 2 - Evaluating the diagrams in Figs. [TT1 
and 1 101 using the expressions for the renormalized diag- 
onal velocities given by Eq. ([5rj]) . one obtains the 
first and second terms in Eq. (|72p respectively. 

This final result completes our comparison showing the 
equivalence between the SBE and the Kubo formalism. 

V. EXCLUDED DIAGRAMS 

In this paper we have considered only diagrams with- 
out crossing disorder paring lines. This approximation is 
self-consistent and justified by Ward identities but incom- 
plete. There is no obvious reason why crossing diagrams 
like the one in Fig. [T^] are small in comparison with oth- 
ers. It is expected, however, that crossing diagrams are 
parametrically different from others and thus the self- 
consistent non-crossing approximation can be justified in 
many cases. 

Another type of terms that we dropped from our dis- 




FIG. 13: A higher order diagram contributing to the skew 
scattering current. 

cussion is due to higher order non-Gaussian correlations. 
An example is shown in Fig. 1131 This diagram leads to 
current that behaves as 1 /rij on the impurity concentra- 
tion. Although formally it is small in comparison with 
similar diagrams from Fig. [11)1 when many different types 
of impurities are involved, the 3rd cumulant of the disor- 
der potential can be suppressed because different types 
of impurities may contribute to it with different signs. 
In contrast, the 4th order cumulants remain of the same 
sign and thus can compete with the 3rd one. We are 
not aware of any work on crossing diagrams and the 4th 
order non-Gaussian skew scattering contributing to the 
AHE. 



VI. CONCLUSION 

We have compared two different techniques commonly 
used for anomalous Hall current calculations. One is 
based on the semiclassical Boltzmann equation (SBE) 
and on the heuristic identification of side-jump and 
anomalous velocity unconventional contributions. The 
second approach is systematic linear response theory, 
which requires evaluation of a Kubo formula. The SBE 
lacks a systematic derivation in a controlled expansion. 
In particular, it is not completely obvious that all un- 
conventional effects have been identified or that these ef- 
fects represent independent contributions that should be 
added together. It is based on the assumption that elec- 
trons in the diffusive regime can be treated as classical 
particles and only the collision term and the expression 
for the velocity need to be calculated quantum mechani- 
cally. This assumption fails, for example, when localiza- 
tion effects are important. The semiclassical approach, 
however, has several advantages, such as the attractive 
transparency of its physical interpretation and its appli- 
cability beyond the linear response approximation. It is 
well known that in the theory of the longitudinal conduc- 
tivity, the SBE is equivalent to the self-consistent non- 
crossing approximation in the Kubo formalism. We have 
extended this correspondence between SBE and Kubo 
approaches by identifying the Feynman diagrams that 
correspond to the unconventional side-jump and anoma- 
lous velocity effects that are important for the anomalous 
Hall effect. We have established a one to one correspon- 
dence between both techniques in leading order of the 
small parameter £ = l/fcpi sc . This direct correspondence 
between particular diagrams and various side-jump and 
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anomalous velocity effects applies only when working in 
the Kubo formalism with the chiral exact Bloch eigen- 
state representation. The identification is based in part 
on the physical transparency of the SBE AHE contribu- 
tions, has been verified by explicit calculation for the 2D 
Dirac model, and is completely general. This connection 
between the two approaches has not been directly es- 
tablished in previous literature, to our knowledge, which 
contains many farraginous results which often appear to 
be inconsistent with each other. 

In addition, in performing this comparison with the 
skew scattering contribution, we have discovered that 
there is a contribution that semiclassically is also related 
to the asymmetric part of the collision term kernel, but 
which leads to a contribution independent on the impu- 
rity concentration and is non-zero even for Gaussian cor- 
relations. The identification of an impurity-independent 
skew scattering contribution has not been made in pre- 
vious literature. Although this contribution to the Hall 
current is independent of impurity concentration, like the 
side-jump and intrinsic contributions, it has a distinct 
physical meaning and thus has generally different depen- 
dence on other parameters, e.g. e F , as one can see from 
our results for the 2D Dirac Hamiltonian. 

The Kubo formalism and the SBE should be consid- 
ered as complementary techniques. The former provides 
a rigorous justification for the latter. The latter provides 
simple semiclassical interpretations of the predictions of 
the former, which is important in developing intuition 
that can be applied to more complicated problems. We 
demonstrated in this work for the first time how very dif- 
ferent approaches of the Hall current calculation converge 
to the same expressions. We hope that the observations 
made in this paper will help increase the consistency of 
various publications devoted to disorder contributions in 
the AHE. 



presence of disorder in the normal spin basis. This ladder 
vertex corrections is defined by v a {e F ) = v a + 5v a (e F ) 



Si a (e F ) = -^^G(e F )J Q G(e F ) (Al) 

k 

+ ^_Y^G R {e F )5J a {e F )G A {e F ) 



We decompose the operators in the Pauli matrix rep- 
resentation Sv a — 8vq<tq + Sv^cii, G — Gq<Jq + GiGi, 
v a = v o °o + v f &ii where the coefficients can be read 
from the expression for G: 

s~<R _ 1 _ eF+iT+v(k x a x +kyi7y) + (A-iri)a z 

U ~ l/G^-S-R ~~ (e F -£ + +ir+)(e F -e-+jr-) 

(A2) 

where T = im t Vg/{Av 2 ), £, = rcos(fl), 7 ± = T ± 
ricos(0), and cos 6 = A/^J (vk) 2 + A 2 . 

We can write a closed equation for the vector 8v after 
multiplying the equation by Ui from the left and taking 
the trace: 



2Sv a = C a + ASv a (A3) 
5v a = (2 - A)~ x C a (A4) 

where the vector C a = (G£ , C£ , C% , C%) and matrix A 
are defined as 
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APPENDIX A: VERTEX CORRECTION 
CALCULATION 

In this appendix we discuss the procedure to calcu- 
late the vertex correction to the velocity operator in the 



C: ^ ^ Yl TrKGJ a G], A ab ^^Y, T ^Ga b G] 

k k 

(A5) 

Here a,b — through 3 and we reserve i,j to label 
x,y,z. Simplifications of these expressions can be ob- 
tained through the use of the following relations 

Tr \&i(Ji> cTj&j'] = Tr[(5u> +iewk(Tk)(Sjf +iejf k '0-' k )] 
= 2(5u<5jj> - SijSfi' + Sij'S^t), (A6) 



Tilaiajcr'j] = Tr[ieij k crkVj'] = 2ie y y. (A7) 



For the case of the Driac model Hamiltonian = v 
and v% = v\ = = 0. Also G x / y cx k x / y and G / z 
do not have any angular dependence. The only non-zero 
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terms in A and G are: 



CI = 



V 



^^Im[G r et Gf v ] = vb 



C y = ^£^ u(G rot G adv_ G ret G adv ) 
k 



va 



V 



A z o — A 0z — 



k 

4mV 2 



V 



]Trc[g 2 g ] 



. _ 2n,V 2 



V 

A-xy Ay X vh 



— ^^(GoGq + G Z G Z — G X G X — G X G X ) 



(A8) 



This gives a block diagonal form to the matrix A so in 
what follows and since Gq = C]f = we will ignore the 
Aoo,A zz ,A z0 ,A() Z terms since irrespective of their value 
Svq = Sv y = 0. The main step we have left to do is to 
compute a and b defined above: 



2mV 2 



V 



^«(G^ ct G§ dv -GTG 



adv \ 



2 ni Vl ^ ^r-> 



6- F + T z - ft 2 - T 2 



V 



D+D- 



(A9) 



r+(4\ 2 + + r 2 _) r_(4A^ + r 2 



x( e | + r 2 -/i 2 -r 2 ), 



(A10) 



where we have defined D± = ((ep — E±) 2 + T 2 t ). Because 
A 2 . = A 2 . = ( 2 F = A and since T± = T(l ± cos 2 9) for 
tp < —h and T± — T(l =p cos 2 9) for tp > h we can write 
this using the fact that 2F/ir = niV^(v+ + v-\. 



a = 4T- 



{e 2 F + T 2 -h 2 -V 2 ) 

(1 + cos 2 6»)(4A 2 +T2sin 2 9)' 

sin 2 9 T 2 

+ 0(- T ) for|e F | > ft (All) 



(1 + cos 



For b we have 
4n,V 2 



where we have used the combination epT = —XT cos 9. 
The matrix A reads (the quadrant of x and y) : 



.1=1 °. b ) andG=( 6 



(A13) 



and putting this back in Eq. IA4l we finally get the renor- 
malized vertex: 

„ (l + cos 2 6<) / sin 2 9 \ 

v y = v8Tcos9 2av2 ax+ \ v + v ii i q Ta\ a v 

A(l + 3 cos 2 9y \ (l + 3cos^f' ' 



(A14) 



APPENDIX B: DERIVATION OF THE OFF 
DIAGONAL CONDUCTIVITY USING CAUSAL 
GREEN'S FUNCTIONS. 

The off-diagonal DC conductivity can be calculated 
using the Kubo formula representation in terms of causal 
Green's functions in the limit to — > 0. We start from 



a xy (w) = — Tr 



de d 2 k 



uj J 2tt(2tt) 2 



r x G{e + u)v y G(e) (Bl) 



with the Green's functions which include a self energy 
related to the scattering from impurities. 

The result of integration over e can be presented in the 
following form 



f(-e)v x 



-2c 



-V 



V (-w)(-2e-w) 



e+ui+v(T-'k+Aa z . -e+vCT-k+Aa 



V 



-2c 



-f{e)vx 



e-\-v(T -k-j-Acfg e—u>+v(T -k+Aer^ 



2£ 



V (-w)(2e-w) 



e+u+v(T-k+Acr z e+v(T-]s.+Aa z 

u)(2e+u) U V 2e 



+f(e)v x - 1 - ; r " ->■ 



+ [f(e + u)-f(e))]u 



\f(e)-f(e + w))] T x 



(2e+w-ir„+ir + ) v v 

t+v(T-k+Aa z -ir +ir + + ir 1 cr z 
2e-ir_+ir + 

c+ui+v(T-k+A(y z +iT a +iT + -iT 1 rj z ^ 



} 



e+v(T-'k+Aa z -iTa+ir + +iT 1 (j z 
2e-ir_+ir + 



b = 



V 



An z V 2 



J2MG r ct Gf v ] 



'kijVj ^ -Th + e F f 



r + (4A 2 + + r 2 _) 



+ 



r_(4A 2 h + r 2 h 



(-r/i + e F r) 



4r CO s6i . r 2 , r , 

TO( — ) for|e F | > ft 



(1 + cos 2 (9)A 



(A12) 



(B2) 

where e = , Note that the renormalization of the ver- 
tex v x — > T x is included only in the last term. Even 
though the 5th term contains the integration near the 
Fermi level, it should include the unrenormalized ver- 
tex. This is because the corresponding Green's functions 
come with different signs of e + lo and e only in the last 
term. For definitcness, we assumed that the Fermi level 
is located in the upper band, tp > A. 

The first four terms in (|B2|1 correspond to the contri- 
bution of states below the Fermi energy, <r™* , which does 
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not have the vertex corrections and is not affected by im- 
purities. The 1st and 2nd terms are associated with the 
filled lower Dirac band, whereas the 3rd and the 4th are 
due to the upper band. In the limit of u) — > it gives 



eVA f d 2 k /(e) - /(-e) 



(2tt) 



(B3) 



The 5th and 6th terms in (|B2[) give two contributions 
to a xy . We can identify them as a h xy and a xy , respec- 
tively, because they are related to integrals of G A G A 
(i.e., with e + lo < 0, e < in a corresponding product of 
causal functions) and G R G A [i.e., with e + uj > 0, e < 
in G(e + u)G(e)]. 

For w, T ->• we find from (|E2|) 



df 



2we W<7 2/ 2e 



~ rLx 2(w+2ir + )e VU V 



t+v(T-V.+Ao z -iTa+r 1 a z - 
2e 



(B4) 



As we see, the second term in (|B4I) gives different results 
in the limits of uj <C T + and a; ^> T + (pure and "dirty" 
cases). In the pure case we should also take Y x — > vcr x . 
Then, as follows from (|B4[) . in the pure limit we get <J xy + 
a xy — 0. It should be noted that each of the terms <r xy 
and a xy is nonzero but they exactly cancel each other. 
For example, calculating a xy we find a xy = e 2 A/4ireF. 
In the limit of oj <C T + , we find from (IF34I) 



^xy'^xy 



e 
4jt 



2 r 



A 



b_Ti_ 



b £o A, 



i - 



(B5) 

where we used the notation = ba x + ca y . After sub- 
stituting the values of b, c, and T we find 



e 2 cos ( 

47T 



1 



8(1 



(1 + 3 cos 2 6) 2 



(B6) 



where cos 6 = A/ejr. Calculating integral (|B3|) we find 

. . nns 

(B7) 



Combining it with (|B5[) we obtain finally 

2e 2 cos 0(1 + cos 2 0) 



■xy 



a xy = 



(B8) 



7T (1 + 3 cos 2 e) 2 
which coincides with Eq.([2]) for V\ = 0. 

APPENDIX C: CALCULATION OF a" y 

In calculating it is useful to consider separately 
integrations over positive and negative e, i.e. a xy = I < + 



/> , where with our simplifications 
-Vx + vr){-^±(G A -{e)- 



-G R -(e)) + i^(G^-(e) - G R ~(e))). 

V'k e k J 



(CI) 



We can get rid of derivative over e by integration by parts 
and noticing that for e < and positive chemical poten- 
tial follows that df(e)/de = 0. Thus 

/< = e 2 def(e) J fa {^f^^S(e - e k ")) 

Further simplifications come from the fact that off- 
diagonal elements of the velocity operator are related 
to interband matrix elements of the coordinate oper- 
ator v+- = (-i[x,H])+~ = iA+-(e+ F - e k ) where 
A + ~ = (tt k \id/dk x \u^) etc. Such elements have the 
property i[A±*A* ± -A± z f : A£ ± ] = F ± , where F+/F- is 
the Berry curvature of the upper/lower sub-band. After 
this 



9 t d 2 k 
1< = ~ e J JW F ~ (k) 



(C3) 



where we used the fact that the lower band is filled and 
hence /(e k ) = 1. 

Calculations of I > is analogous, 

I>=-€l +Q °def(e)J^(vt-v y +- 
-^x + ^-)(^-tA gA+ V G R+ (z))+ (C4) 

k 

-_i (G^+(e)-G^(e))) 

v e k E k I 



however for positive e follows that df(e)/de = —S(e — 
ep) 7^ and thus extra terms will appear after integration 
by parts. 



/> = /; + 1» (C5) 



where 



r> = -e 2 



'eZ<e F 



d 2 k 

WT 2 



F+(k) 



(C6) 



and 



> 4-tt J— oo ae J (2-^)^ 



-:- ^-^<- ( ^ +(£) _ G A +(£))] 



(C7) 



Comparing (|C7[) and (fSU)) we find that they differ only 
by sign thus we can write that 
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We conclude that a^ y contains a term that exactly can- 
jj 2 . , A eels the intrinsic term ai% nt) from a 1 



(C8) 

TP-n,\\ _I(mt) 
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